Crossover between Mean-Field and Ising Critical Behavior in a Lattice-Gas 

Reaction-Diffusion Model 
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Lattice-gas models for CO oxidation can exhibit a discontinuous nonequilibrium transition be- 
tween reactive and inactive states, which disappears above a critical CO-desorption rate. Using 
finite-size-scaling analysis, we demonstrate a crossover from Ising to mean-field behavior at the crit- 
ical point, with increasing surface mobility of adsorbed CO or with decreasing system size. This 
behavior is elucidated by analogy with that of equilibrium Ising-type systems with long-range in- 
teractions. 



I. INTRODUCTION 



Nonequilibrium phase transitions and pattern forma- 
tion occur in a broad variety of physical, chemical, biolog- 
ical, and sociological systems>i*^ Such complex behavior 
sometimes resembles more familiar phenomena in equilib- 
rium systems. However, there are also some fundamental 
differences due in part to the lack of free energy minimiza- 
tion or detailed-balance constraints for nonequilibrium 
systems. The latter feature also provides a challenge 
for development of a rigorous theoretical framework. 
Nonetheless, the spatial-temporal behavior of nonequi- 
librium transitions can be effectively studied using one of 
two paradigms. The first is through a set of partial dif- 
ferential equations (e.g., reaction-diffusion equation for 
chemical reactions), where phase transitions are related 
to bifurcations of the underlying nonlinear dynamicsii^ 
By default, criticality is governed by mean-field (MF) be- 
havior, although Ginzburg criteria can be developed to 
characterize noise-dominated regimes^ The second ap- 
proach utilizes microscopic interacting particle systems 
of which nonequilibrium lattice-gas models form an im- 
portant subclass. For these models, one can successfully 
apply techniques such as finite-size-scaling, and concepts 
including that of universality classes familiar from critical 
phenomena in equilibrium phase transitions." 

Many types of nonequilibrium phase transitions have 
been studied: continuous transitions to absorbing states^ 
(which are often in the universality class of directed per- 
colation), discontinuous transitions between reactive and 
inactive states,^'® and order-disorder transitions,^ all of 
which occur in adsorption-desorption or reaction models; 
driven lattice-gases^"® (where the driving breaks detailed 
balance); dynamic Ising models with both spin fiip and 
spin exchange dynamics, but at two different tempera- 
tures leading to interesting crossover phenomena.^ For 



some of these systems, it is possible to derive rigorous 
reaction-diffusion equations in the hydrodynamic limit 
of rapid exchange or diffusion. ^'^^ 

111 this paper, we consider specifically surface reaction 
models, which incorporate adsorption-desorption and re- 
action steps, as well as surface mobility. Our motivation 
is that these models not only contain the essential ingre- 
dients to elucidate complex phenomena observed in real 
surface catalysis, but also exhibit many fascinating fea- 
tures of the above idealized models. On the one hand, 
they exhibit complex dynamics and phase transitions as- 
sociated with the reaction mechanism (and perhaps with 
adspecies interactions), but also surface diffusion rates 
provide a key parameter which can tune behavior. In 
fact, while many theoretical studies have considered the 
regime of limited or zero mobility, it is the hydrodynamic 
regime of rapid mobility that is most often realized in 
experiments.^ 

Our focus here is on first-order phase transitions and 
associated critical phenomena in models for CO oxida- 
tion on metal surfaces, behavior analogous to that of 
the ferromagnetic Ising model. The central question is 
"how does critical behavior depend on surface mobility" , 
anticipating that MF behavior could well apply in the 
hydrodynamic regime, but that a crossover may occur 
to non-MF behavior for lower mobility. Through precise 
finite-size-scaling (FSS) studies of simulation data, we 
confirm this prediction and quantify crossover behavior. 
Wc also comment on experimental realization of crossover 
for CO oxidation in high-pressure nanoscale systems. 



II. REACTION MODEL: SPECIFICATION AND 
BEHAVIOR 

We now describe our lattice-gas reaction-diffusion 
model, where gas (ads) denote gas-phase (adsorbed- 
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phase) species. The following steps are implemented 
using a square lattice of adsorption sites with periodic 
boundary conditions. (i) CO(gas) adsorbs on single 
empty sites at rate p = pco (chosen between and 
1) and desorbs with rate d. (ii) 02(gas) adsorbs dis- 
sociatively onto a pair of diagonally adjacent empty sites 
at rate P02 = 1 ~ PcOi provided all six neighbors are 
free of O(ads). This "eight-site" rule reflects strong 
nearest-neighbor (NN) 0(ads)-0(ads) repulsions. Since 
O(ads) is treated as immobile, this adsorption rule en- 
sures that O(ads) never occupy adjacent sites, (iii) ad- 
jacent CO (ads) and O(ads) react at rate k (set to unity 
here), (iv) CO(ads) hops to adjacent empty sites at rate 
h. T his model has also been discussed elsewhere, e.g., 
Refs.m 

Conventional kinetic Monte Carlo (KMC) simulations 
are used to assess model behavior for finite h. Noting 
that mobility of CO(ads) is often very high under ultra- 
high vacuum conditions, we also perform a direct anal- 
ysis of limiting behavior for h = 00 using a "hybrid" 
treatment: here the distribution of O(ads) is described 
within a lattice-gas framework, but one only tracks the 
number of CO (ads) and assumes that they are randomly 
distributed on sites not occupied by 0(ads)i& For very 
large h and low 0-coverages, this is valid. 

We now briefly review the steady state behavior of this 
model for an infinite system. First, consider the case of 
finite h < 00.^'^ For low d, one typically finds a first-order 
transition at some p = pco = P* between a reactive state 
with low CO-coverage (6*00) (for p < p*) and an inactive 
state with high (^co) (for P > P*)- The transition dis- 
appears as d increases above some critical value dc{h). 
Second, consider the hybrid model with h — 00. Here 
one finds a region of bistability with both reactive and 
inactive states (the discontinuous transition at p = p* for 
finite h corresponding to the equistability point). Bista- 
bility disappears at a cusp bifurcation upon increasing 
d above some dc{oo). A coherent picture for both cases 
of infinite and finite h comes from the observation that 
decreasing h decreases the degree of metastability or hys- 
teresis in the system. 

Next, we describe how behavior changes for finite L. 
For the reaction model with finite h, the discontinuous 
transition in the (6*00) versus pco mentioned above is 
rounded for finite L, but becomes sharper as L increases. 
The trend is analogous to behavior for the equilibrium 
Ising model in finite systems. In the hybrid model with 
h — 00, there is also a smooth transition in the (Oco) 
versus pco for finite L. This refiects the feature that 
for L < 00, the system can make noise-induced jumps 
between the low {6co) reactive and high {6co) inac- 
tive states, and that the relative weight of these states 
changes smoothly with pco ■ This transition also becomes 
sharper with increasing L, occurring at the equistability 
point for L — > 00. See Ref. [h] for details. 

Finally, to place our study in a broader context, we 
note that our CO oxidation model corresponds to a mod- 
ified version of the Ziff-Gulari-Barshad (ZGB) model, ^ 



and its extension with desorption of CO (ads). The key 
modifications in our model are inclusion of: a) hopping of 
CO (ads); b) NN exclusion of O(ads); and c) finite rather 
than infinite reaction rate. All these features are impor- 
tant for realistic modeling of CO oxidation. However, it 
is primarily the first feature [hopping of CO (ads)] which 
impacts the critical behavior of the reactive-inactive tran- 
sition studied in this paper. Furthermore, it is clear 
that our conclusions about critical behavior would ap- 
ply to the ZGB model, modified to include desorption 
and hopping of CO (ads). The second feature [NN exclu- 
sion of O(ads)] results in an oxygen poisoning transition 
in the ZBG model^ being replaced by an order-disorder 
transitionii This does not affect critical behavior of the 
reactive-inactive transition, with one caveat. For small 
h, NN exclusion does also lead to a loss of the reactive- 
inactive transition (see Sec. II V|) . The third modification 
to the original ZGB model is that instead of fc = 00, 
we choose k = 1, i.e., the reaction rate equals the total 
adsorption rate. The choice is of course somewhat arbi- 
trary, but the basic behavior of these models does not 
change varying the reaction rate from 0(1) to 00. For 
further discussion of effects of var ying k on the steady 
state bifurcation diagram, see Ref. Il2i 



III. CRITICAL POINT DETERMINATION AND 
FSS ANALYSIS 

Unlike the Ising model, transitions in reaction model 
do not involve simple symmetry-breaking. As in the 
liquid-vapor phase separation problem, one needs to lo- 
cate the critical point in a two parameter space, i.e., (p, d) 
for the reaction model. However here one also has the 
disadvantage that computationally efficient techniques 
for equilibrium systems (e.g., histogram-reweighting and 
cluster algorithms) do not apply. Thus, when using nu- 
merical techniques such as FSS, careful analysis of sim- 
ulation data is necessary. Below, we briefiy describe our 
procedure. 

There are several reasonable ways to define the effec- 
tive transition point p ~ p*^ for finite L. One natural 
choice is the point where the change in the relevant or- 
der parameter (i.e., the CO coverage {Oco)) is greatest. 
All choices should converge to the same value as L di- 
verges (for fixed d < dc), and should yield the same 
critical exponents in the FSS analysis described below. 
For a given d < d^ we measure CO-coverage (^co) ver- 
sus p for two different system sizes, say L and 2L. We 
find that a convenient definition of is the point where 
the curves cross, i.e., where {0co{p*l))l = (^'co(pI))2L- 
The technique is similar in its underlying motivation to 
a method for equilibrium first-order transitions by Borgs 
and Jankeii^ 

After determining the transition pressure p* for each 
d value, one can study the critical behavior using the 
usual FSS techniques. For example, the quantity xl = 
-^^((^co) ~ (^co)^), which is related to susceptibility in 
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equilibrium systems, is assumed to have the following 
behavior 
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(1) 



if the transition belongs to the Ising universality class. 
Here v = 1 and 7 = 7/4 are the critical exponents for the 
correlation length and susceptibility respectively (in 2D). 
In contrast, for the MF universality class, one ha^i^ii^ 
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where Z? = 2 is the spatial dimension of the system so 
that is the volume (or area in our case) of the system. 
Defining Rl = X2l/xl, m either case, one has 



RL{dc) ^2''', as L — * oo 



(3) 



where the size scaling exponent 7 is -D/2 = 1 for the MF 
universality class, and 7/1^ = 7/4 for the Ising universal- 
ity class. 

Since one does not know the value of dc a priori, a 
convenient way to determine both the critical point and 
exponent is by finding the crossing point of Rhid) and 
R2L{d), so that 



Rl (dc ) = R2L {di 



271. 



(4) 



7l can be considered as the effective critical exponent for 
finite systems (see Sec. lIV|l . Note that in determining d^ 
and 71,, we need simulations of system of linear sizes L, 
2L, and 4L. 

Figure^shows Rl versus d for the hybrid model. The 
FSS argument above predicts that Rl for different L's 
will cross at {dc, Rc) where Rc = 2 for MF criticality, and 
= i^l^ « 3.364 for Ising (in 2D) criticahty. Figure [H 
clearly shows MF behavior for system up to L = 128. 
Note that in the limit of L ^ 00, i?^ is a step function 
with Rl = 4 when d < dc, and Rl = 1 when d > dc.— 

Table |2 lists the effective critical exponent obtained 
using the above procedure for the hybrid model. All es- 
timates for 7i are close to unity, consistent with the pre- 
diction of MF universality class where 7 = Extrap- 
olation to i = cxD assuming a 1/i^ finite size correction 
(consistent with MF criticality) gives dc — 0.05258(5) 
and Pc — 0.41327(5). Assuming a 1/L finite size cor- 
rection for 7i gives 7 — 0.995(6). It is possible that 
corrections to scaling are described by other exponents. 
However, our analysis shows that the range of system 
sizes is large enough so that the value of dc and the con- 
sistency with MF behavior is not dependent on the form 
of the corrections. Note that MF behavior occurs despite 
the presence of spatial correlations in the distribution of 
O(ads) due to limited mobility of O (ads) and interactions 
between 0(ads)i 

Figure|21shows Rl versus d for the reaction model with 
finite h = 1. The crossing point appears to approach 2^/** 
as L increases, consistent with the prediction of Ising uni- 
versality. Limitations in analysis of larger system sizes 



2.20 



2.10 



2.00 r 




0.055 



FIG. 1: Ratio, Rl, of mean-square fluctuation amplitudes 
for systems of size 2L and L at equistability points for the 
hybrid model. L ranges from 8, 16, 32 to 64 with increasing 
steepness. The dotted line plot Rl = 2, the MF prediction. 



TABLE I: Effective critical point and critical exponent for 
fluctuations of ^co for the hybrid model obtained from finite- 
size-scaling analysis. 



L 


dc 




JL 


8 


0.05304(6) 


0.4138(1) 


0.948(2) 


16 


0.05275(5) 


0.41343(8) 


0.968(4) 


32 


0.05264(3) 


0.41330(6) 


0.981(5) 




0.024 



FIG. 2: Ratio, Rl, of mean-square fluctuations amplitudes 
for systems of size 2L and L at equistability points for the 
reaction model with flnite CO diffusion model h= 1. L ranges 
from 8, 16, and 32 with increasing steepness. The dotted line 
shows Rl = 2^"''*, the Ising prediction. 



precludes definitive convergence of FSS results for either 
the hybrid model (Fig. ^ or finite h (Fig. [21 However, 
comparing these cases reveals very different behavior, 
and we argue that the above assignment of universality 
classes is quite reasonable and natural. 
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TABLE II: Variation of the critical desorption rate dc{h) with 
CO hop rate h. Results are obtained from FSS scahng using 
system sizes L —16, 32, and 64. 



h 


0.5 


1 


2 


4 


10 


20 oo 


dc{h) 


0.012 


0.019 


0.027 


0.035 


0.042 


0.047 0.0527 



From the slope of Rl at the crossing point and the 
scaling form x[{d — dc)L^] in Eqs. and Q one can 
also estimate the exponent related to the rounding of the 
critical region due to finite sizes. Using results in Fig. 
we find that 9 — 1.0(1) for the largest system sizes, con- 
sistent with the MF value 9 — D /2 — 1. Uncertainties 
in the data of Fig. |5] are too large to obtain accurate es- 
timate of 9, but they are consistent with the Ising value 
l/v = 1. In two dimensions, 9 can not be used to distin- 
guish between the MF and Ising universality class. 

We note one previous study by Tome and Dickmanii 
of the ZGB lattice-gas model^ for CO oxidation (O2 ad- 
sorption on adjacent sites, infinitely fast reaction, no CO 
diffusion) modified to include CO-desorption. They mea- 
sure the shift of the critical point with system size and 
found that dc{L) — dc{oo) ^ L^^ where A = 1 consis- 
tent with results of the two-dimensional Ising model. In- 
deed, using a similar method, we found that A ~ 1 for 
h — 1 and A w 2 for the hybrid model. This result 
is consistent with our conclusion regarding universality 
classes. However, care must be taken with the interpre- 
tation of the shift exponent. A, because it is sensitive 
to the boundary conditions (which were chosen to be 
periodic in both the above studies) in the case of MF 
universality where hyperscaling is violated. For exam- 
ple, choosing free boundary conditions (which modifies 
adsorption and reaction processes near the boundaries) 
could lead to modified A for the hybrid model. 



IV. CROSSOVER BEHAVIOR 

Table nil shows the variation with h of the critical point 
dc for the CO poisoning transition obtained using FSS. 
The h = 00 value is taken from the hybrid model. By 
analogy with equilibrium studies, discussed below and in 
Ref. 14, we assume that the shift of dc{h) away from the 
limit of ft. ^ cx) scales as l/h (although there could be 
logarithmic corrections). 

As an interesting aside, we note that extrapolation be- 
havior to the regime of small h suggests that dc{h) van- 
ishes ai h = ht ~ 0.2. Thus, for small Q < h < ht, the 
CO poisoning transition is continuous for d = 0, and does 
not exist for d > Q. This is in contrast to the ZGB model 
where a first-order CO poisoning transition exists even 
for immobile CO. For the case of d = and h < ht, our 
epidemic analysis (cf. Ref. 0) indicates that the tran- 
sition to the absorbing (CO poisoned) state, belongs to 
the directed percolation universality class. 




FIG. 3: Effective critical exponent 7 versus h. See text for 
details. 



Thus, to summarize our results for the critical point 
in our reaction model with fc = 1, we find that the crit- 
ical desorption rate dc increases from dc ~ for (small) 
h — ht to dc — 0.0527 for h = 00. In contrast, for the 
ZGB model (no NN exclusion of O) with k — 00, modi- 
fied to include CO desorption and CO mobility, one finds 
that dc = 0.0 1^^1^^ ioi h — increasing to dc — 2/3 for 
h — ooi^ Further refining this ZGB-type model to in- 
corporate k = 1 (rather than k — 00), one finds that 
dc = O.OaSa for ft = increasing to dc = 0.142 for 

Of more central interest to this paper is the variation 
with h of the (effective) critical exponents. Upon in- 
creasing the CO diffusion rate, ft, how does one crossover 
from Ising criticality (applying for finite h > ht) to the 
MF criticality of the hybrid model (which corresponds to 
first taking the limit of ft — > 00). A related question is: 
for finite (possibly large) ft, how does critical behavior 
depend on finite system size? 

In the literature on crossover studies, effective expo- 
nents are usually defined as the slope of the log-log plot 
of the measured critical quantities versus deviation from 
the critical point, e.g., jeS = — dln^/dln \d — dc\. A dif- 
ferent definition, which is more amenable to numerical 
simulations is to fix the parameters at the critical point, 
while changing the system size L. The effective (reduced) 
exponent can then be defined as 7cff = dlnx/dlnL. In 
principle, this derivative (or its finite difference approxi- 
mation) should be evaluated at the critical point for the 
infinite system, dc- In practice, dc is unknown a pri- 
ori, and estimation by extrapolation produces additional 
errors (particularly using simulation data for limited sys- 
tem size). Thus instead, one just uses dc obtained from 
the FSS analysis in Eq. Q, and therefore also from 
Eq. I0J as the effective exponent. 

Figure |2Ia) shows 7^ versus variable 1/ft for ft = 0.5, 
1, 2, 4, and 10. The system size is indicated in the fig- 
ure. As Fig. I^Ia) shows, 7 approaches 1 as ft —> 00, and 
approaches 7/4 when ft is small. This behavior can be un- 
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derstood in terms of the theoretical framework presented 
below. 

Phenomenological analysis of nonequilibrium dynam- 
ics (and more specifically of critical phenomena) in 
Hamiltonian systems with nonconserved order parame- 
ters has traditionally been formulated in terms of time- 
dependent Ginzburg-Landau equations (TDGLE) M- For 
multi-component order parameter, 0q,, the TDGLE has 
the general form 



^a(x,t) 

dt 



= -r- 



SF 



Ca(x,i), 



(5) 



where ^q, is Gaussian white noise terms. Here F is 
the coarse-grained free energy functional, which for an 
Ising-type model (with single-component order parame- 
ter) with interactions of range R can be written asi^ti^ 



F= d: 



(6) 



The quartic potential V{(j)) converts from double- to 
single-well form with increasing temperature through 
the critical point. This type of formulation extends to 
nonequilibrium reaction systems characterized by a vec- 
tor c of adspecies concentrations where the TDGLE are 
replaced by reaction-diffusion equations (RDE) of the 
form 

9c(x, t)/dt cx -dVes{c)/dc + V • D • Vc + C- (7) 

Analogous to V, the effective potential Vcs converts from 
double- to single-well passing through the critical point 
via the bistable region. In our problem, the diffusion 
tensor satisfies D oc h (the microscopic hop rate) except 
for small hm^ Thus, comparing the RDE with TDGLE 
after performing the functional differentiation, it is clear 
that h^^^ plays the role of R. This result is reasonable 
if one recognizes these types of reaction-diffusion mod- 
els (with finite reaction rate and typically large diffu- 
sion rates) exhibit a characteristic diffusion length, Ldiff, 
which scales like Ldiff ~ {h/Ky^"^. Here, h is the micro- 
scopic hop rate, and K is some effective rate for the over- 
all adsorption-reaction process. See Refs. 1^123. Thus, 
refining the above statement, we can say that Ldifi plays 
the role of R. Finally, we also note that from analysis of 
either the TDGLE or the RDE for infinite systems, it is 
long known that fluctuations dominate close enough to 
the critical point for spatial dimension D < A. 

Recently, extensive studies have been performed on fi- 
nite size effects in equilibrium Ising-type systems with 
long-range interactions. From scaling and renormaliza- 
tion procedureSfi^ii^ as well as numerical simulations j^k 
it has been shown that the crossover from MF to Ising 
behavior in finite L x L systems is governed by the ratio 
L/R?. Thus, by analogy, one expects that the crossover 
for reaction-diffusion systems is governed by L//i (in 2D). 
However, equilibrium and reactions-diffusion systems are 
quite different in detail: the Hamiltonian formulation. 



which is central to analysis of the former, does not ex- 
ist for the latter; the microscopic nature of diffusion is 
distinct from that of long-range interactions. Thus, it is 
appropriate to test numerically the above scaling hypoth- 
esis. In Fig. |2Ib), we replot the data in Fig. |2Ia) using 
scaling variable L/h. Data collapse, while imperfect for 
L = 8, is quite good for L > 16, supporting our proposal 
that that L/h is the correct size scaling variable. 

Of course, it would be preferable to set the phenomeno- 
logical analysis on a more rigorous footing. In general, 
the strategy of extending mean-field-typc formulations to 
include noise via stochastic partial differential equations 
can be unreliable. An alternative rigorous approach is to 
attempt to map the exact master equations for the model 
onto a stochastic field theory.^^ This strategy has been 
applied for a variety of toy reaction-diffusion models, 
e.g., demonstrating that the introduction of diffusion can 
modify the number of adsorbing states and the universal- 
ity class in pair-contact models with diffusion.^^ However, 
these models are quite different from our reaction-limited 
monomer-dimer type models, where the utility of the rig- 
orous approach has yet to be demonstrated. 



V. CONCLUSIONS 

Catalytic CO oxidation on noble metal surfaces pro- 
vides a convenient (and important) chemical system for 
which lattice-gas modeling can be used to rigorously as- 
sess criticality in a reaction-diffusion system. We have 
systematically studied the role that CO diffusion plays 
on the criticality of the CO poisoning transition. We find 
a crossover from MF criticality for fast diffusion [which 
occurs despite spatial correlations in the distribution of 
O(ads)], to Ising criticality for limited diffusion. The nat- 
ural variable describing this crossover is proposed based 
on phenomenological arguments. The proposition is sup- 
ported by numerical simulations. 

Recent experiments on CO oxidation in nanoscale sys- 
tems can probe fluctuations and critical behavior iSSiSS 
Thus, the considerations of this paper, including the 
study of crossover as a function of system size, become 
particularly relevant. Furthermore, the crossover behav- 
ior described in this study may be realized physically 
as one makes the transition from low-pressure to high- 
pressure catalysis. The latter could produce a decrease 
by many orders of magnitude of the diffusion rates (rel- 
ative to other rates) that results from the higher surface 



coverages 
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